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In order to investigate the mechanism of ferromagnetic ordering in the new n-type magnetic 
semiconductor (In,Fe)As co-doped with Be, we have performed X-ray absorption spectroscopy and 
X-ray magnetic circular dichroism (XMCD) studies of ferromagnetic and paramagnetic samples. 

The spectral line shapes suggest that the ferromagnetism is intrinsic originating from Fe atoms 
incorporated into the Zinc-blende-type InAs lattice. The magnetization curves of Fe measured by 
XMCD were well reproduced by the superposition of a Langevin function representing superparam- 
agnetic (SPM) behavior of nano-scale ferromagnetic domains and a T-linear function representing 
Curie-Weiss paramagnetism even much above the Curie temperatures. The data at 20 K showed 
a deviation from the Langevin behavior, suggesting a gradual establishment of macroscopic ferro¬ 
magnetism on lowering temperature. The existence of nano-scale ferromagnetic domains indicated 
by the SPM behavior suggests spatial fluctuations of Fe concentration on the nano-scale. 
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I. INTRODUCTION 


a 


Since the discovery of ferromagnetism in (In,Mn)As [l| and (Ga,Mn)As [2|, the Mn-doped III-V ferromagnetic 
semiconductors (FMS) have been intensively studied, driven by the scientific interest of how ferromagnetism is 
realized by such a small amount of magnetic impurities and by future technological applications utilizing 

both the charge and spin degrees of freedom SHTj. Although their electronic structures and the mechanism of the 


ferromagnetism still remain highly controversial 8l-ll3l|. it is widely believed that carriers mediate ferromagnetic 


interaction between spatially separated magnetic ions. This kind of ferromagnetism is called carrier-induced ferro¬ 
magnetism, and its microscopic elucidation and utilization have been the main subjects in the field of semiconductor 
spintronics. In the case of Mn-doped III-V compounds, however, there is a limitation that only p-type materials 
are available because the Mn^'*' ions substituting the group-III elements work not only as local magnetic moments 
but also as acceptors providing the system with hole carriers. 

Recently, Hai et al. 


n 


14l-ll8j| have synthesized a new n-type III-V FMS, namely, Fe- and Be-codoped InAs. Usually, 


the substitution of Fe^"*" ions for the In sites does not provide charge carriers, resulting in Curie-Weiss paramagnetism 
19|. However, Hai et al. have shown that by co-doping (In,Fe)As with Be atoms, which enter the interstitial sites 
and act as double donors, electrons are introduced into the system and mediate ferromagnetic interaction between Fe 
atoms. Here, unlike the Mn-doped materials, (In,Fe)As:Be has the advantage that one can control the concentration 
of magnetic ions and that of charge carriers independently. 

In order to reveal the nature of the ferromagnetism of this new FMS, it is necessary to characterize the electronic 
and magnetic structure on the microscopic level. For this purpose, we have performed X-ray absorption spec¬ 
troscopy (XAS) and X-ray magnetic circular dichroism (XMCD) measurements. XMCD is defined as the difference 
between the absorption of right- and left-handed circularly polarized X-rays and is a powerful method to probe the 
element-specific electronic and magnetic properties of compounds by utilizing the absorption edges of constituent 
elements. Furthermore, one can exclude extrinsic contributions to magnetism such as magnetic contaminations and 
the diamagnetic response of the substrate in the case of thin hlms and that of the host material in the case of 
FMS. Here, such diamagnetic contributions prohibit the extraction of the intrinsic paramagnetic compon ent 
which is important to characterize the complex magnetism of FMS. In addition, XMCD sum rules 
possible to obtain the orbital and spin magnetic moments separately. A recent XMCD study 23| have shown that 
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221 make it 


the orbital moment of (In,Fe)As is signihcantly larger than that of Fe metal. Since Fe^"*" (d®) does not have orbital 
magnetic moment, the observed finite orbital magnetic moment would represent the charge-transfer d^L states of 
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, where L denotes a ligand hole. 


In the previous study 15|, ferromagnetic domains of '^10 /rm were observed at sufficiently low temperatures well 
below Curie temperature (Tc) by magneto-optical imaging, explaining the hysteretic behavior of magnetization. 
At high temperatures much above Tc, however, magnetization still strongly depends on magnetic held, which 
implies the existence of rather small ferromagnetic domains of ^nm acting as a superparamagnet. Under such 
circumstances, it is important to know how the ferromagnetism evolves as a function of temperature and magnetic 
held on a more microscopic level, i.e., on the nano-meter scale. In this paper, we have performed a careful analysis 
of XMCD data in order to address the issue of the evolution of magnetism from high to low temperatures. 


II. EXPERIMENT 


Three samples Ino.gsFeo.osAsiBe (sample 1), Ino.gFeo.iAsiBe (sample 2) and Ino. 95 Feo. 05 As (sample 3) without 
Be were synthesized using the low-temperature molecular beam epitaxy (LT-MBE) method as follows. The 20-30 
nm-thick (In,Fe)As layers were grown on InAs(OOl) substrates at 240°C after growing 50 nm-thick InAs buffer 
layers at 500° C. In order to prevent surface oxidation, the samples were covered by sub-nm-thick amorphous 
As capping layers. The amount of Be was estimated to be 2.6 x 10^® cm“®. Samples 1 and 2 were found to 
be ferromagnetic, whose Curie temperatures were 15 K and 40 K, respectively, determined by the Arrott plot 
of visible magnetic circular dichroism (vis-MCD) intensities, while sample 3 was paramagnetic. Careful sample 
characterizations by transmission electron microscopy (TEM), energy dispersive X-ray spectroscopy (EDX), and 
three dimensional atom probe (3DAP) proved that there were no Fe related precipitates or secondary phases in the 
samples 1^. Experiments were done at the undulator beamline BL-16A2 of Photon Eactory (PF), High Energy 


Accelerator Research Organization (KEK). Temperature was varied from 20 K to 260 K and magnetic field from 0 
T to 5 T. The total electron-yield (TEY) mode was used. The direction of the incident X-rays and magnetic field 
was perpendicular to the sample surfaces, and XMCD spectra were obtained by changing the polarization of X-rays 
while the direction of magnetic field was fixed. Backgrounds of the spectra were assumed to be the summation of 
a linear function and two-step inverse tangent functions representing the Fe Tg.a edge jumps 241. 


III. RESULTS 


The XAS and XMCD spectra of (In,Fe)As with and without Be are shown in Fig. [T] together with the spectra 
of Fe metal (ferromagnetic) and y-FegOs (ferrimagnetic) as references. Clear XMCD signals were obtained for the 
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FIG. 1. XAS (a) and XMCD (b) spectra of (In,Fe)As in comparison with those of metallic Fe [2^ and 7 -Fe 203 [1^. The 
XAS and XMCD spectra have been normalized so that the peak intensities of XAS are equal to one. For the (In,Fe)As 
sample with 10% Fe, the spectra have been decomposed into the intrinsic and extrinsic components as depicted with red 
dashed and blue dashed curves, respectively, (c) Fe Ls-edge XAS and XMCD spectra of the 10% Fe-doped (In,Fe)As sample 
at different temperatures and magnetic fields. The XMCD spectra as well as XAS spectra have been normalized to its peak 
height. In order to emphasize weak structures, the second derivative of the XAS spectra with reversed sign are also shown 
by gray solid curve. 


ferromagnetic (In,Fe)As samples and virtually no signal was detected for the paramagnetic (In,Fe)As sample as 
expected. The XAS spectra at the Fe A 3 edges are mainly composed of two structures a and /3. The small shoulder 
structure (3 could be attributed to an extrinsic origin of Fe oxides at the surface [23| , considering that its energy is 


close to the peak energy of 7 -Fe 203 and that its intensity changes in non-systematic ways with varying magnetic 
field and temperature, as shown in Fig. [TJc). This probably resulted from the change of the beam position between 
different scans and inhomogeneous distribution of Fe oxides on the surface. Utilizing this non-systematic spectral 
change, we have deduced intrinsic XAS and XMCD spectra in the following manner 26|. First, using two spectra 


Si and S 2 with different intensities of structure /3, the intrinsic spectrum Sint was obtained as Sint cx Si — pS 2 , where 
p was chosen so that structure (3 vanished. Then, the extrinsic spectrum Sext was obtained as Sext oc Si — gSint, 


where q was chosen so that line shape of Sext agreed with that of 7 -Fe 203 . The decomposed intrinsic and extrinsic 
spectra for the 10% Fe-doped sample are shown in Figs. [TJa) and[TJb). It was difficult to apply the same procedure 
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to the 5% Fe-doped sample because structure /3 is less prominent and the intensity of structure /3 did not change 
significantly between different scans. However, using the intrinsic and extrinsic spectra for the 10% sample, we have 
estimated the mixing ratio of the intrinsic and extrinsic components, that is, how much oxides exist at the surface. 
Thus we could obtain the correct values of the magnetic moments for the 5% Fe-doped sample as well as 10% 
Fe-doped one using XMCD sum rules [2l|,[23. We note that, except for the extrinsic shoulder structures /3, the line 
shapes of the XAS spectra were independent of the Fe and Be contents and shared a broad single-peak structure. 
This means that the local electronic structures of Fe atoms is basically the same between the three samples, notably 
between the ferromagnetic and paramagnetic ones. The result that multiplet structures reminiscent of the localized 


28, 3 


3d systems [27| are absent suggests that the d electrons of Fe are relatively delocalized [28|, l29j probably through 
the strong hybridization between the Fe 3d orbitals and the ligand s,p orbitals. 

The deduced intrinsic XAS and XMCD spectra of (ln,Fe)As:Be resemble those of Fe metal at a first glance, 
and therefore, one might suspect that the ferromagnetism originates from Fe metal precipitates. However, this 
possibility can be ruled out from following reasons. First, the XAS spectra of the paramagnetic sample showed the 
same spectral line shape as those of the ferromagnetic samples. If the observed spectral line shapes of (In,Fe)As were 
due to Fe metal particles, the XMCD intensity should have been comparable among the three (In,Fe)As samples 
including the paramagnetic one. Second, looking at the magnified spectra of the Fe L 3 edge presented in Fig. [He), 
the peak positions of the XAS and XMCD do not coincide with each other unlike the XAS and XMCD spectra of 


Fe metal 


30|. Corresponding to the peak in the XMCD spectra, there is a tiny shoulder a' in the XAS spectra 


0.5 eV below its peak position a. This shoulder is clearly seen in the second derivative of the XAS spectra with 
respect to photon energy, but less prominent in the spectra of Fe metal 30|. From these observations, we conclude 


that the ferromagnetism in this system is distinct from that of Fe metal and of the intrinsic origin arising from Fe 
atoms incorporated into the Zinc-blende-type InAs lattice. 

In order to obtain the magnetic moment of Fe for the ferromagnetic samples, we have applied the XMCD sum 
rules 21 , 221 : 


4 f. Act dw 

^ 2,3 


^^2,3 

^orb — f j '^h'i 

3 r cr acj 

^ ^ 2,3 

(1) 

6 C Act dw — 4 f. Act dw 


1 w i-'3 ^-‘-'2,3 

^spin H” t'fTlT — r j 

L (7 dcj 

^ ^ 2,3 

(2) 

CT = CT+ -I- CT_, Act = CT+ — ct_ , 

( 3 ) 


where morb and rrispin are the orbital and spin magnetic moments in units of ps /atom, respectively, and ttit is the 
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TABLE I. Spin and orbital magnetic moments of Fe in (In,Fe)As:Be and Fe metal. 



‘^orb/^spin 

TTZorb 

^spin 

Iiio.QsFeo.osAsiBe 

0.088 ± 0.016 

0.22* 

2.49* 

Ino.gFeo.iAsiBe 

0.075 ± 0.015 

0.21* 

2.80* 

Ino.gsFeo.osAsiBe [23] 

0.065 ± 0.014 

0.10 

1.60 

Fe bcc [2^ 

0.043 ± 0.001 

0.085 

1.98 


*measured at fJ,oH= 5 T and T = 20 K 


expectation value of the magnetic dipole operator which is negligibly small for an atomic site with high symmetry 
such as Td or Oh 3l|. cr+ and (t_ denote absorption cross sections for X-rays with positive and negative helicity, 


respectively, Uh the number of 3d holes. We applied these equations to the decomposed intrinsic spectra instead of 
the raw spectra. The area of the intrinsic XAS spectra, namely the denominator of the equations, were estimated 
to be 20% smaller than that of the raw XAS spectra for the 10% Fe-doped sample, and ^ 5% smaller for the 
5% Fe-doped one. As for the XMCD spectra, the extrinsic components were almost absent and did not contribute 
to the area of the XMCD spectra and hence to the numerator of Eqs. (1) and (2). Note that the existence of the 
extrinsic components in the XAS spectra only changes the calculated spin moment mgpin and orbital moment morb 
by the same factor while keeping the ratio morb/wspin unchanged. 

Table I summarizes the orbital and spin moments of (In,Fe)As:Be deduced using the XMCD sum rules and those 


of bcc Fe metal 


24| . For the calculation of the moments, we assumed the number of electrons in the Fe 3d states 


32l | . Being consistent with the previous 


to be 5 and the correction factor for the spin sum rule to be 0.685 for Fe^’*' 
report 23|, the ratio of the orbital moment to the spin moment niorb/wspin, which is not affected by the assumption 


of n/i = 5 and the existence of the extrinsic component in the XAS spectra, is larger than that of Fe metal and is 
positive. This can be explained if the dPL configuration is mixed into the ground state due to charge transfer from 
ligand orbitals to the Fe 3d orbitals. 

Magnetization (M) per Fe atom at various magnetic fields and temperatures has thus been deduced from the 
XMCD intensities and are plotted in Figs. Ha) and Hb) ■ Here, the horizontal axis is the effective magnetic field 
corrected for the demagnetizing field of the thin films although this correction is small in the case of the diluted 


magnetic systems. In the present experimental set-up, hysteresis, which was observed in the previous study 23|, 
was too small to be detected. It is worth noting that the M-H curves were almost identical between 5% and 10% 
Fe-doped samples even though the Tc’s were different. At low magnetic fields, a steep increase of magnetization was 
observed even above Tc and, at high magnetic fields, linear increase of magnetization was observed even at the low 
temperature of 20 K where the ferromagnetic component would be saturated. This behavior can be understood if two 
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FIG. 2. XMCD intensities of (In,Fe)As:Be samples with 5% and 10% Fe doping, (a), (b) Magnetic field dependence of 
magnetization deduced from the Fe L 2 , 3 -edge XMCD intensities at various temperatures, (c), (d) The same data as (a) 
and (b) plotted against fioH/T. At the low temperature of 20 K, the magnetization shows a clear deviation from the high 
temperature data. In panel (a)-(d), solid curves represent the fitting curves by the summation of the superparamagnetic 
Langevin function and the paramagnetic linear function. The paramagnetic components are separately shown by dashed 
lines, (e), (f) Slope of the magnetization curves at zero magnetic field against ^qH/T, namely AM/A |jr->o = x'l' s-s 

a function of temperature. 


contributions to the magnetization are superimposed, namely, superparamagnetism from nano-scale ferromagnetic 
domains (FM/SPM domains) and paramagnetism from isolated magnetic moments. Therefore, we have fitted the 
following summation of the Langevin function L(^) and a T-linear function to the experimental XMCD intensity: 


C T f CfloH 

L(^) = coth (0 - 


(4) 

(5) 


where x {0 < x < 1) denotes the fraction of Fe atoms participating in the superparamagnetism or ferromagnetism, 
fj, the total magnetic moment per nano-scale FM/SPM domain, C the Curie constant of the Fe^’*' ion with 5 fiB, 
and Ta the Weiss temperature representing the antiferromagenetic interaction between paramagnetic Fe^"*" ions. 


Figures [2Dc) and[2Dd) show the magnetization curves at various temperatures replotted as functions of iiqH/T. 
They approximately fall onto a single curve, which is a fingerprint of superparamagnetism, except for the data at 
20 K close to Tq- The deviation may represent the onset of macroscopic ferromagnetism at 20 K. Figures [2De) and 
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FIG. 3. (a) Fraction x of Fe atoms participating in ferromagnetism or superparamagnetism to all the Fe atoms, (b) Total 
magnetic moment /r inside one FM/SPM domain, obtained by fitting Eq. (4) to the data, (c) Density of FM/SPM domain 
N estimated from the fitting parameters x and fj,. 


[Iff) show the slopes of magnetization at zero magnetic field with respect to ^qH/T, i.e. AM/A = X^- 

If the system follows the Langevin law of superparamagneitsm, yT should be constant. Hence the drop of xT 
around 50 K on decreasing temperature indicates the gradual establishment of macroscopic ferromagnetism from 
the nano-scale superparamagnetism. 


The M-H data were fitted by Eq. (4) assuming Ta = 30 K, which is nearly equal to the Weiss temperature of 


paramagnetic (Ga,Fe)As without carrier doping (= 32 K) 


19|, and the results of the fit are shown in Figs. [2ja)-[2jd) 


by solid curves. In the figures, the paramagnetic components are separately shown by dashed lines. Note that a finite 
value of Fa larger than 30 K was necessary to fit to the data. Figures[3l^a) and[3jb) show the temperature dependence 
of fitting parameters x and /i. The fraction of Fe atoms participating in the ferromagnetism or superparamagnetism 
(x) is estimated to be 40%-60% and increasing as the temperature is decreased. The magnetic moment per FM/SPM 
domain (fi) are 200-300 ^b, which corresponds to 40-60 Fe atoms, and they decrease as the temperature is decreased. 
Interestingly, the estimated values are very similar between the two samples. The estimated density of FM/SPM 
domains JV are shown in Fig. He) and found to be ~ 10^® cm comparable to the density of Be atoms. 
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FIG. 4. Schematic picture of the formation of nano-scale FM domains in (In,Fe)As:Be. There exist FM Fe-rich domains and 
paramagnetic isolated Fe ions. Here the density of Fe atoms and the size and the number of FM domains reflect the obtained 
value for the 10% Fe-doped sample. 


IV. DISCUSSION 


There may be two possible scenarios to explain the formation of nano-scale FM/SPM domains observed for 
(In,Fe)As. One is the bound magnetic polaron (BMP) model 33|, [3^, where relatively localized carriers stabilize 
ferromagnetic spin alignment extending over the scale of nm. In this model, the polarons overlap with each other and 
form macroscopic ferromagnetic domains at the critical concentration where the system satisfies the Mott criterion 
for insulator-to-metal transition. Second is the formation of Fe rich regions during the growth while maintaining 
the zinc-blende crystal structure, often discussed in the context of spinodal decomposition [3^, l36| , as schematically 


shown in Fig. SI Considering the metallic behavior of (In,Fe)As 15|, in particular, the small effective mass of 


3 


electrons in the conduction band [la, the BMP model is less likely because it predicts a heavy effective mass for 
the charge carriers, and the second scenario is more probable. As noted above, the M-H behaviors are almost 
identical between the two ferromagnetic samples with different Fe contents, meaning one FM/SPM domain consists 
of 40-60 Fe atoms for both samples. If the first scenario is correct and the distribution of Fe atoms is rather uniform, 
electrons are bound to charged impurities such as Be^“'" to form magnetic polarons and there should be difference in 
the total magnetic moment of a FM/SPM domain between the two samples because the number of Fe atoms within 
the same Bohr radius of electrons should be different. On the other hand, in the case of the second scenario, where 
Fe-rich FM/SPM domains are formed in the InAs host, the observed independence of magnetic moment on the Fe 
content can be explained if Fe-rich regions with optimized number of Fe atoms are formed. Suppose an extreme 
case where Fe atoms occupy all the In sites inside the FM/SPM regions, the radius of the domain containing 40-60 
Fe atoms is estimated to be ^ I nm, which is difficult to detect by conventional x-ray diffraction. From the fact 
that the density of FM/SPM domains are comparable to that of Be atoms, one can speculate that the FM Fe-rich 
domains were preferentially formed around Be impurities. Actuall y, i t was reported that the co-doping of shallow 


impurities affects the distribution or aggregation of magnetic ions 


37|. 
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We shall discuss below how the magnetism in this system behaves with temperature. As can be seen from Figs. 
ISJa) and[3Jb), the lower the temperature becomes, the more Fe atoms participate in ferromagnetism but the smaller 
the average size of FM/SPM domains becomes. This behavior can be understood if there exists size distribution 
of Fe rich regions. As described in Fig. |4l at high temperatures, only large Fe-rich regions are ferromagnetic. On 
lowering temperature, smaller Fe-rich regions becomes ferromagnetic as described in the middle panel of Fig. 01 and 
hence the average magnetic moment of FM/SPM domains (/x) decreases [Fig. [H[b)]. On the other hand, the fraction 
of Fe atoms participating in ferromagnetism or superparamagnetism (x) increases [Fig. |31[a)] and the density of 
FM/SPM domains (N) increases [Fig. Eljc)]. Further decreasing the temperature below Tq would result in the 
overlap of the FM/SPM domains as shown in the left panel of Fig. 01 and the macroscopic ferromagnetism, which 
cannot be described by the Langevin law, starts to be stabilized. The reason for the drop of xT at 20 K [Figs. [21[e) 
and[2l[f)] instead of a divergence with the emergence of macroscopic ferromagnetism may be that the ferromagnetism 
is soft in (In,Fe)As. For example, let us write the magnetization of a system where superparamagnetism and soft 
ferromagnetism coexist as 


= (1 - ) + yta.nh{H/Ho), 


M,. 


ksT 


( 6 ) 


where the first term and second term represent the superparamagnetism and the soft ferromagnetism, respectively. 
y denotes the fraction of ferromagnetism and Hq is the saturation field. The second term is just an empirical 
I, but one can assume that the magnetization is independent of temperature as long as the temperature 


expression 


is sufficiently low. Using this expression, yT can be written as 




(7) 


Therefore, if yofi/SkB is larger than T/Hq, i.e. Hq > 3kBT/yoy,, yT would show a drop with the emergence of 
ferromagnetic term. 

At the present stage, one can conclude that the large magnetization of the present (In,Fe)As:Be samples above Tq 
originates from the nano-scale FM domains. The previously observed ferromagnetic domain of ~10 ym at sufficiently 
low temperature would contain a large number of 10^ assuming the domain size of lO^m x lO^m x 20nm) the 
nano-scale FM domains revealed in the present work, and would result from ferromagnetic coupling between the 
FM/SPM domains and/or magnetic anisotropy of the FM/SPM domains. 
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V. CONCLUSION 

We have performed XMCD measurements on the new n-type ferromagnetic semiconductors (In,Fe)As:Be with 
different Fe and Be contents. Despite the fact that Tq was different between the samples, the XAS and XMCD line 
shapes and magnetization data were almost identical. The analysis of the spectral line shapes and the magnetic 
moment using sum rules suggest that the ferromagnetism is intrinsic originating from the Fe atoms incorporated 
into the Zinc-blende-type InAs lattice, and not from segregated secondary phases. From the magnetization curves 
at various temperatures even much above Tq, a SPM-like steep increase in the low H region and a paramagnetic 
linear increase in the high H region were observed, which implies the co-existence of FM domains of ^nm size 
and paramagnetic isolated Fe^^ ions. The fraction of Fe atoms participating in the ferromagnetism or superpara¬ 
magnetism was estimated to be 40-60% and the magnetic moment per FM/SPM domain to be 200-300 ^b- The 
formation of the FM/SPM nano-domains can be explained by the formation of Fe-rich regions during the material 
growth possibly around Be atoms. 
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